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Abstract. A vector field a on a Riemannian manifold M is said to be har- 
monic if there exists a member of a 2-parameter family of generalised Cheeger- 
Gromoll metrics on TM with respect to which a is a harmonic section. If M 
is a simply-connected non-flat space form other than the 2-sphere, examples 
are obtained of conformal vector fields that are harmonic. In particular, the 
harmonic Killing fields and conformal gradient fields are classified, a loop of 
non-congruent harmonic conformal fields on the hyperbolic plane constructed, 
and the 2-dimensional classification achieved for conformal fields. A classifica- 
tion is then given of all harmonic quadratic gradient fields on spheres. 



1. Introduction 

Let a be a vector field on an n-dimensional Riemannian manifold (M, g). It was 
observed in [9] and [8] that if M is compact and a is a harmonic map [6] into TM 
equipped with the Sasaki metric h [12] then a is parallel; furthermore this remains 
the case under the less stringent and arguably more natural condition that a is 
a harmonic section of TM; ie. a critical point of (vertical) energy with respect 
to variations through vector fields [TS]. This rigidity may be overcome for vector 
fields of unit length by further restricting the class of variations to unit vector 
fields; in other words, the study of harmonic sections of the unit tangent bundle. 
These harmonic unit Gelds have acquired an extensive literature; see for example 
the bibliography of [TJ. Unfortunately its scope is restricted by the topology of M. 
The essential reason for the rigidity of the energy functional in this context is that 
the restriction of the Sasaki metric to tangent spaces is flat, creating an analogy 
between harmonic sections of TM and harmonic (vector-valued) functions on M; 
whereas harmonic sections of the unit tangent bundle are analogous to harmonic 
maps into spheres [13]. To complete the analogy, ordinary derivatives are replaced 
by covariant derivatives, and the Laplace-Beltrami operator for functions by the 
rough Laplacian for vector fields; thus the equations for a harmonic unit field are: 

V*V(T = |Vcr| 2 cr, (1.1) 

where: 

V*Vct = -trV 2 a = -EiV^^CT, (1.2) 
for any local g-orthonormal tangent frame {Ei} on M. 
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In [3] it was proposed to address the rigidity problem for vector fields by em- 
bedding h in a 2-parameter family h p>q (p, g 6 K) of generalised Cheeger-Gromoll 
metrics, within which h = /io,cb the Cheeger-Gromoll metric [5] appears as 
and /i2,o is the stereographic metric. The h Pjq all belong to the infinite-dimensional 
family of g-natural metrics on TM [TJ [2], but are much more tightly controlled, 
being constructed from a spherically symmetric family of metrics on K™ via the 
Kaluza-Klein procedure. Then u is said to be (p, q)-harmonic if er is a harmonic 
section of TM with respect to h Ptq (the metric g on M is fixed throughout). From 
[3] the corresponding Euler-Lagrange equations are: 

T p , q (a) = T p (a) - (f> p . q (a)(7 = 0, (1.3) 

where: 

T p (a) = (1 + \a\ 2 )V*Va + 2pV VF a, 
<f> p , q (a) = p|W| 2 - pq\VF\ 2 - g(l + \a\ 2 )AF. 

Here F is an abbreviation for \ \cr\ 2 , VT denotes the gradient vector field, and AF 
the Laplacian. If a is a unit field then it is easily seen that (|1.3[) reduces to 
with p = 2; thus harmonic unit fields are (2, q)-harmonic, for all q. If <r is parallel 
then a is (p, g)-harmonic for all (p, q). 

A notable feature of the generalised Cheeger-Gromoll family is that if q < then 
h p>q has variable signature. Precisely, if for r > the ball and sphere subbundles 
of TM are denoted: 

BM(r) = {Ie TM : \X\ < r}, SM(r) ={Ie TM : |X| = r], 

and g < 0, then h Piq is Riemannian on the q-Riemannian ball bundle BM(r) where 
r = y/—l/q, degenerate (in radial vertical directions) on SM(r), and Lorentz- 
ian on the exterior of BM(r). However, for all r > y/—l/q the sphere bundle 
SM(r) is a spacelike hypersurface of (TM, h Ptq ), and the induced Riemannian met- 
ric is vertically homothetic to that induced by the Sasaki metric. We say that a 
is q-Riemannian if the pullback cr*h Ptq is Riemannian. This condition, which is 
independent of p, is equivalent to either a having constant length, or: 

q\o(x)\ 2 > -1, forallxeM. (1.5) 

Note that (|1.5j) allows a to touch the boundary of the g-Riemannian ball bundle. 

In this paper we say that a is a harmonic vector field if a is (p, q)-harmonic for 
some pair (p,q), which we refer to as metric parameters for a. Under the weaker 
condition that T p (a) is pointwise collinear with a we say that a is p-preharmonic. If 
the metric parameters for a harmonic a are unique we will say that er is metrically 
unique. The previously mentioned examples show that this need not be the case. 
However, with one exception, all the examples considered in this paper turn out to 
be metrically unique; and the exception has just two pairs of metric parameters. 
Nevertheless, it is interesting to learn that metric non-uniqueness can occur for 
vector fields of non-constant length. Part of the problem when looking for harmonic 
vector fields is to identify metric parameters, and to this end the following result 
of [3] provides some guidance. 
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Proposition 1.1. Suppose M is compact and a is a harmonic vector field of non- 
constant length. For each p there exists at most one pair of metric parameters (p, q) 
for a. Furthermore if \p\ ^ 1, or p > 1 and ||cr||oo ^ W-P ~ 1> then q < 0, with 
q<l-p/2 lfp^2. 

Proposition 1 1.1 1 does not assert the existence of a lower bound for q on compact 
manifolds, or that a harmonic vector field is necessarily g-Riemannian. However this 
turns out to be the case for all the examples considered in this paper. Interestingly, 
the contrapositive can be used to provide bigger than expected lower bounds for 
their image diameter. 

A feature of the harmonicity equation (11.31) that is implicit in Proposition II . II is 
its non-invariance under dilations of TM . In fact, if a is harmonic and r € R then 
unless r = 0, — 1, or the length of a is constant, it cannot be inferred that ra is har- 
monic, even under a change of metric parameters. Our examples illustrate this very 
clearly. Thus harmonic vector fields come with a preferred scale, and determining 
this is another element of the problem. It also means that the space of harmonic 
vector fields is not self-evidently contractible. However, equation (|1.3[) is invariant 
under the action of the isometry group of (M, g). Thus if a is (p, q , )-harmonic then 
so is any vector field congruent to a, and harmonic vector fields therefore need only 
be classified up to congruence. (This could justifiably be considered a basic require- 
ment of the theory) . Continuous families of harmonic vector fields may therefore be 
easily constructed by the action of 1-parameter subgroups of isometrics. However in 
this paper we exhibit a 1-parameter family of non-congruent harmonic vector fields. 
Such families appear in general when (M,g, J) is a Kahler manifold, in which case 
there is a circle action: 

e lt .a — (cosi)cr + (sint) Ja, 

with respect to which T Ptq is equivariant: 

T Ptq (e lt .(j) = e lt .T VA (a). 

Thus, if a is (p, g)-harmonic then so is e %t .a for all t. Borrowing the terminology 
of classical minimal surface theory, we refer to this S^-orbit of a as the associate 
family of harmonic vector fields, and Ja as the conjugate harmonic field. 

Throughout the paper Af is a simply-connected non-flat real space form; ie. the 
sphere S n or real hyperbolic space H n . Our approach is unashamedly extrinsic, 
which we find provides greater geometric insight. As far as possible we try to treat 
the two cases simultaneously, although there are often marked differences between 
them, and certain situations where the hyperbolic case requires more attention. 
Underlying our search for harmonic vector fields are a number of unifying themes. 
One is that all vector fields considered can be constructed from conformal gradient 
Gelds. Such fields were considered as the prime example in [3J, but for the sphere 
only, and it was shown that whilst none exist on the 2-sphere, for each dimension 
greater than two there exists a unique (up to congruence) harmonic representative, 
which is also metrically unique. In Section [2] we extend that classification to hy- 
perbolic space. Notable differences are the appearance of a 2-dimensional example, 
and a family of metrically non- unique examples in higher dimensions. Many of the 
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basic computational elements set out in Section [2] will be re- used in our analyses of 
more complicated vector fields. In Section |3] we classify all harmonic Killing fields 
on M . It came as somewhat of a surprise to us to learn that Killing fields can be 
constructed from conformal gradient fields, and we utilise such decompositions later 
in the paper. However in Section [3j whilst decomposition into conformal gradients 
affords an attractive analysis of harmonic Killing fields on spheres, in the hyperbolic 
case there are additional complications that make the procedure rather unwieldy; 
we have therefore chosen an alternative approach that works well for both. We first 
show that all harmonic Killing fields on hyperbolic space are infinitesimal rotations, 
thus eliminating infinitesimal translations and fields of parabolic type. These, along 
with all spherical Killing fields, are classified up to congruence by their spectrum 
of rotational frequencies, which we refer to as the field's twists, and we then show 
that harmonicity demands that these are all equal, or in other words that the field 
is balanced. Finally, after discounting the Hopf fields on odd-dimensional spheres 
(which are the only Killing fields of constant length, and well-known to be har- 
monic |14| 116]). we show that there is a unique value for the twist, and determine 
the metric parameters, which are also unique. In the spherical case it transpires 
that harmonic Killing fields of non-constant length exist only in dimension four 
and higher. However harmonic Killing fields exist on hyperbolic spaces of all di- 
mensions; in particular, on the hyperbolic plane, we note the appearance of the 
conjugate of the unique (up to congruence) harmonic conformal gradient field. In 
Section U we examine loxodromic vector fields. These conformal fields provide a 
natural deformation between conformal gradient fields and certain types of Killing 
fields: non-Hopf fields on the sphere, or infinitesimal rotations on hyperbolic space. 
On the hyperbolic plane, the associate family of the harmonic vector fields uncov- 
ered in Sections [2] and [3] is loxodromic. We show that these are in fact the only 
harmonic loxodromic fields in any dimension. In Section[5]we study another natural 
conformal deformation between conformal gradients and Killing fields, that is ap- 
plicable to infinitesimal hyperbolic translations. This deformation passes through 
a dipole Geld, which is of interest in its own right. However, the only possible har- 
monic vector fields in the entire deformation are the "endpoints" ; therefore in the 
hyperbolic case only the conformal gradient field "endpoint" can be harmonic, and 
on the 2-sphere there are no harmonic vector fields of this type. In Section [6] we 
consolidate results of Sections [4] and [5] by classifying all harmonic conformal vector 
fields in 2-dimensions, confirming non-existence on the 2-sphere. Finally in Section 
[3 with the 2-sphere somewhat in mind, we classify all harmonic quadratic gradient 
fields on spheres. In fact, these exist only on odd-dimensional spheres of dimension 
five and higher. Each such sphere has a unique (up to congruence) harmonic repre- 
sentative <7, which is also metrically unique. As one might hope, the geometry of a 
is pleasingly simple: a pair of equidimensional mutually orthogonal totally geodesic 
submanifolds of maximal dimension play the role of source and sink, and a flows 
along the inter-connecting great circles with a specified maximum speed. 

A second unifying theme of the paper is that all our examples of harmonic vector 
fields turn out to be eigenfunctions of the rough Laplacian. In this case a is clearly 
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O-preharmonic, and p-preharmonic for p ^ if and only if: 

V VF a = Co, (1.6) 
for some smooth function ( : M — > R, where necessarily: 

H 2 c-|v^| 2 . (1.7) 

Equation (jl.6|l may be regarded as a partial integrability condition for the full 
harmonicity equation (|1.3p . and is notable for its independence of the metric pa- 
rameters. We say an eigenfunction a is preharmonic if (|1.6j) holds, and refer to C 
as the spinnaker of a. Use of (|1.7|) together with the Weitzenbbck identity: 

(V*Vo,cr) = |Vaf + AF (1.8) 

then reduces the harmonicity equation (|1.3[> to the following PDE for F: 

(p + q + 2qF)AF + 2p(l + qF)C + v(l + 2(1 - p)F) = 0, (1.9) 

where v is the eigenvalue. 

2. CONFORMAL GRADIENT FIELDS 

Throughout the paper we use the hyperboloid model of hyperbolic space. Thus, 
let R™' 1 denote equipped with the Lorentzian inner product: 

(x,y) = xiyi H h x n y n - x n+x y n+1 , for all x,y e R n+1 . 

Then: 

H n = {xe M™' 1 : (x,x) = -1, x„+i > 0}. 

For the sphere we use the standard model: 

S n = {xe « n+1 : (x,x) = 1}, 

where (-,-) now denotes the Euclidean inner product on R™ +1 . In both cases the 
Riemannian metric g on M is that induced by (-,-), and for notational convenience 
we will often also denote g by (-,-). In order to handle spherical and hyperbolic 
computations simultaneously we introduce the indicator symbol e, which is set to 
1 when M — S n and —1 when M = H n . Furthermore, we let V denote ambient 
(n + l)-dimensional Euclidean or Lorentzian space, as appropriate. 

Let a G V be any vector, and set p = (a, a). If V is Lorentzian then p may 
be negative, in which case we assume without loss of generality that a is future 
oriented. Let a: M — > R be the restriction of the covector metrically dual to a: 

a{x) = (a, x), for all x e M. (2.1) 

The conformal gradient Geld A on M with pole a is then defined a = Va, the 
spherical/hyperbolic gradient of a. We sometimes refer to a as the potential of o. 
It follows from the linearity of a that for all X g TM: 

(a, X) = da(X) = a(X) = (a, X) . (2.2) 

Therefore if {Ei} is an orthonormal basis of T X M: 

(r{x)=^2 i (a(x),E i )E i =^2 i (a,E i )E i = a-ea{x)x, (2.3) 
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since {E\, . . . , E n , x} is an orthonormal basis of V. In particular: 

\<t\ 2 = (i-ea 2 . (2.4) 

In the spherical case a has precisely two zeros (viz. ±a/ \ffl), and /x is the maximum 
value of \a\ 2 , which is attained on the equatorial hypersphere orthogonal to a. In 
the hyperbolic case it follows from (|2.3|) and (|2.4|) that a has unbounded length. 
Furthermore a has a zero only when a is timelike, in which case the integral curves of 
a are, as point sets, the geodesic rays from the zero (viz. a/y/—fi). If a is spacelike 
then the integral curves of a are the ultraparallel family of geodesies orthogonal 
to the totally geodesic hypersurface a 1 ' D H n , which we refer to as the equator for 
ex. Furthermore fj, is the minimum value of \cr\ 2 , which is attained on the equator. 
If a is lightlike the integral curves of a are asymptotically parallel geodesies, and 
a — > asymptotically. In the spherical case, and the hyperbolic cases with fj, ^ 0, 
a is determined up to congruence by [i. 

The covariant derivative of M is characterised by the Gauss formula [10] : 

V x Y = D x Y + e(X,Y)x, (2.5) 

where D denotes the standard directional derivative of M n+1 , and x is interpreted 
as a unit normal to M. Hence by (|2.3|) and (j2.2[) : 

V x cr = — ea(X)x — ea(x)X + e(X, a)x = —eaX. (2.6) 

From ()2.6j) the divergence of a is: 

diva = Ei^a.Si) = -ena. (2.7) 

Furthermore the second covariant derivative is: 

from which, by (II. 2| : 

V*V(T = £(7. 

So cr is an eigenfunction of the rough Laplacian. 

Lemma 2.1. All conformal gradient fields on a non-flat space form are prehar- 
monic, with spinnaker £ = s(/i — 2F). 

Proof. From (gH) : 

VF = -m\7a = -eacr, (2.9) 

hence by ([276]) : 

V VF a = -eaVF = a 2 er. 
Therefore cr is preharmonic, and ^ follows by comparison with (|2.4[) . □ 

Theorem 2.2. Lei cr be a non-trivial conformal gradient field on the non-flat space 
form M, with pole a € V and (a, a) = /!. 

(1) If j-i (ie. M = S n , or M = H n and a has no zeros) then a is harmonic 
if and only if n > 2 and [i = l/(n — 2). Furthermore a is metrically unique, with 
metric parameters (p, (?) = (n + 1, 2 — n). 



= -e( ( r,X)y, 

(2.8) 
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(2) If fi < (ie. M = H n and a has a zero) then a is harmonic if and only if 
fi = — 1. Furthermore if n = 2 then a is metrically unique, with metric parameters 
(3, —1/2), and if n > 2 then a has precisely two pairs of metric parameters: 

(a) (p,q) = (n + 1,1 -n + l/n); (b) (p, q) = (1/(2 - n), 0). 

Proof. By (2l|), ([HI) and (12^41) : 

AF = -div VF = ediv(acr) = e\a\ 2 - no? = s{2{n+ 1)F - np), 

where we have used the general Riemannian identity: 

div(fX) = fdWX + df(X). (2.10) 

Using this and Lemma [2. 11 the harmonicity equation (|1.9|) may therefore be written: 

(p + q + 2qF){2{n + 1)F - np,) + 2p{l + qF)(fi - 2F) + 2(1 - p)F + 1 = 0. 

This is polynomial in F, the coefficients of which yield the following three equations: 

(n + l-p)? = 0, (2.11) 
(n-2)p+(n + l)q+(p-n)q(i = -l, (2.12) 
((2-n)p-nq)fi= -1. (2.13) 

If q = then (2 — n)p = 1 and p, = — 1. Otherwise p = n + 1, whence adding (|2 . 12[) 
and p,13p yields: 

(n-2 + 9 )(l + /i) =0. 

If p = — 1 then n<7 = 1 + n — n 2 , by (|2.13p . Otherwise q = 2 — n, and (n~2)p = 1 
by (j2~T3)) . □ 

In view of Proposition ll.il the metric parameters (1/(2 — n), 0) cannot occur on 
the n-sphere. All other harmonic conformal gradient fields a have metric parameter 
q < 0. As remarked in [3], in the spherical case cr is g-Riemannian, with cr(M) 
touching the boundary of the Riemannian ball bundle on the equator of a. In 
the hyperbolic case the opposite occurs: if a has no zeros then a(M) lies entirely 
outside the Riemannian ball bundle, and touches its boundary on the equator of a. 
The case n — 2 is distinctive, for by [U Theorem 2.8 and Theorem 3.5] the metric 
parameters (3, —1/2) endow the fibres (resp. total space) of the Riemannian ball 
bundle with positive sectional (resp. scalar) curvature. 

3. Killing Fields 

Let cr be a Killing field on M . Then a is an eigenfunction of the rough Laplacian; 
for by a general curvature identity for Killing fields: 

V*Vcr = Ric(cr) = s{n- l)a. (3.1) 

For computational purposes we view a as the restriction to M of a skew-symmetric 
linear transformation A of V; thus: 

a(x) = A{x), for all x € M. (3.2) 
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In the spherical case, the following result may be deduced without computation 
from the Divergence Theorem. However the hyperbolic case requires the Laplacian 
AF, which in both cases will be required for later results. 

Proposition 3.1. Every harmonic Killing field a on M is preharmonic. 

Proof. We need to show that, if non-trivial, a cannot have metric parameter p = 0. 
It follows from (11.31) . (jl.4[) and (|3.1[) that a is (0, q)-harmonic if and only if: 

qAF = e(l - n). (3.3) 
To compute AF we first apply the Gauss formula (|2.5j) to (13. 2ft : 

V x cr = A(X) + s{X, a)x = A(X) - e(A(X),x)x. (3.4) 
Then, summing over an orthonormal tangent frame {Ei}: 
VF = (V E a,a)Ei = (A{E t ) , A{x)) E t 

= -A 2 {x) + e{A 2 {x), x)x = -A 2 (x) ~ e\a\ 2 x. (3.5) 
Using (|2.5[) once again: 

V X (VF) = D X (VF) +e(X,VF)x 

= -A 2 (X) + e(A 2 (X),x)x - e\a\ 2 X, (3.6) 

from which: 

AF = - div VF = -(V^VF),^) = en\a\ 2 - (A(F 4 ), A(Ei)) 

= £ (n + l)|a| 2 -L4| 2 , (3.7) 

where \A\ 2 denotes the Lorentzian norm when e = — 1 (sec (|3.11l) below; however it 
suffices to note that this is constant). If a is (0, g)-harmonic, comparison of Q3.3P 
with (|3.7[) shows that a has constant length, which contradicts (|3.3j) . □ 

Before proceeding, we pause to summarise some relevant geometric features of 
Killing fields on space forms. In the spherical case there exists an orthonormal 
basis {ei} of Euclidean R™ +1 , a positive integer r ^ (n + l)/2 and real numbers 
uji ^ ui2 ^ ■ ■ ■ ^ uj r > such that for all 1 ^ i ^ r and 2r < j ^ n + 1: 

A(e 2 i~i) = oJie 2 i, A(e 2 i) = -u>ie 2 i-i, A(ej) = 0. (3.8) 

The string (wi,...,w r ) determines A (as a linear endomorphism) , hence a (as a 
vector field), up to congruence. Geometrically, the flow of a is a superposition of 
rotations around the r orthogonal great circles cut out by the 2-planes e 2 i-i A e 2 i 
of M n+1 , with angular frequencies wi, . . . , ui r - We refer to r as the rotational rank 
of <7, and the Ui as the twists of a. In particular, uj\ = Halloo- If 2r ^ n then 
kerA intersects S n in a great (n — 2r)-sphere, which we refer to as the axis of a. 
If u)\ = • • • = L) r = u) then we say that a is balanced with twist lj. In this case a is 
congruent to wS r where: 



Y, r (x) = (-x 2 ,x x , . . . ,-x 2r ,x 2r -x,Q, ...,0). 



(3.9) 
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If n is odd and 2r = n+ 1 then E r is the standard Hopf field; otherwise S r restricts 
to a Hopf field on the great (2r — l)-sphere orthogonal to its axis. 

In the hyperbolic case, fix w £ H n and set v = A(w) = <r(w) £ T w H n . Let 
R = R w : V — > V be the orthogonal projection of A into T w H n , defined: 

R(u) = A(u + (u, w)w) + (A(u + (u, w)w), w)w 
= A(u) + (it, w)v — (u, v)w, 

for all u £ V. Then R(w) — and the restriction of R to T w H n is a skew-symmetric 
endomorphism of a Euclidean vector space, which may therefore be put into normal 
form 032]), with 2r ^ n. Now define T = T w = A - R w : 

T{u) — (u,v)w — (u,w)v. (3.10) 

Since R and T are skew-symmetric transformations of V, their restrictions to H n 
are also Killing vector fields, which we refer to as the rotational (resp. translational) 
part of o at w. We also refer to r as the rotational rank of a at w, and u>x,...,oj r 
as the twists of a at w. The dependence of the rotational rank and twists on w is 
constrained as follows. Let denote the indefinite inner product induced on 

the space of Lorentzian skew-symmetric endomorphisms of V: 

(i4i,A 2 )i=tr(Aio(A 2 )t), 

where A^ denotes the Lorentz adjoint. In terms of a Lorentz-orthonormal frame 
containing w: 

(A 1 ,A 2 ) L - Ei(Mei),A 2 (ei)} - (A^w), A 2 (w)), (3.11) 
and it follows that: 

\A\l =\R\l -\T\l^2r^-2r\ 

where r = \v\ and ui is the quadratic mean twist: u 2 — (uj 2 + • ■ ■ + ui 2 )/r. (The 
second equation follows from the invariance of (-, -) j, under Lorentzian congruences.) 
Thus the function rw 2 — r 2 is constant over H n . Its sign determines the nature of 
a. If positive, then a is an infinitesimal rotation. In this case ker(A) is timelike, 
and ker(A) n H n is a totally geodesic submanifold of H n which we refer to as the 
axis of rotation. The translational part of a vanishes at all points on the axis. 
Furthermore the rotational rank and twists of a are invariant along the axis, and 
determine a up to congruence. In particular, if oji = ■ ■ ■ = oj r = oj along the axis 
then we say that a is balanced with twist uj and rotational rank r, in which case 
<7 is congruent to ojS r with S r defined by (|3.9I) . However since 2r < n + 1, in 
contrast to the spherical case S r has at least one zero, and has precisely one zero 
when 2r = n. We refer to E r , in both the spherical and hyperbolic 
generalised Hopf field. If ru 2 ^ r 2 then ker(A) is spacelike or lightlike, so a has no 
zeros. When the inequality is strict, a is an infinitesimal translation. In this case 
there is a unique geodesic of H n that is an integral curve of cr, which we refer to as 
the axis of translation. The rotational part of a vanishes at all points on the axis, 
and the length of a along the axis, which is necessarily constant, determines a up 
to congruence. Finally, if ruS 2 ~ r 2 then a is of parabolic type. 
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Our next step, in the light of Proposition 13. 1[ is to explore the implications of 
preharmonicity. We therefore compute, using ()3.4j) and ()3.5j) : 

V VF cr = A(VF) - e(A{VF), x)x = -A 3 {x) - e\a\ 2 a. (3.12) 

(Note that since A 3 is also skew-symmetric, its restriction to M is indeed a (Killing) 
vector field on M.) Thus a is preharmonic if and only if there exists a smooth 
function A: M —> K such that: 

A 3 (x) - X(x)A(x) = 0, for all x eM. (3.13) 

Differentiation of (|3.13l) yields: 

A 3 {X) - X(x)A(X) = dX(X)A(x), for all X e T X M. (3.14) 

It follows from (|3.13[) and (|3.14|) that for each x G M the linear map A 3 — X(x)A 
has rank at most one. But non-trivial skew-symmetric transformations of V have 
rank at least two; so A 3 — X(x)A = for all x € M. Since A{x) vanishes at worst 
on the intersection of M with an (n — l)-dimensional subspace of R™ +1 — which 
has measure zero — it then follows from (|3.14[) that dX = generically, hence by 
continuity identically. Therefore by connectedness A is constant. We conclude that 
a is preharmonic if and only if: 

A 3 - XA = 0, for some A € M. (3.15) 

Comparison with (|3.12l) shows that the spinnaker is: 

C = -(A + £ |a| 2 ). (3.16) 

We first explore the implications of this for the sphere. 

Proposition 3.2. Let a be a Killing field on S n . Then a is preharmonic if and 
only if a is balanced. 

Proof. It follows from (|3. 151) that a is preharmonic if and only if im(A) is an eigen- 
subspace of the symmetric operator A 2 . Since A is skew-symmetric with respect to 
a Euclidean metric, im(^4) and ker(A) are (orthogonal) complementary subspaces 
of V. In particular A, and hence A 2 , restricts to an isomorphism of im(A), which 
is therefore the direct sum of all the non-zero eigenspaces of A 2 . Thus a is prehar- 
monic if and only if A 2 has a unique non-zero eigenvalue. From the normal form 
(|3.8p the non-zero eigenvalues of A 2 are — lu 2 , 1 ^ i ^ r. Hence a is preharmonic if 
and only if a is balanced. □ 

Corollary 3.3. Let a be a Killing field of maximal rotational rank on an odd- 
dimensional sphere. Then a is preharmonic if and only if a is a constant multiple 
of a Hopf field. Therefore if a is preharmonic then a is harmonic. 

We now scrutinise hyperbolic space. Choose w G H n and decompose A = R + T , 
to obtain the rotational and translational parts of a at w. Set v — A(w) € T w H n . 
Our next result is diametrically opposite to the conclusion of Corollary 13. 31 

Proposition 3.4. Let a be an infinitesimal translation of H n . Then a is prehar- 
monic, but not harmonic. 
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Proof. Suppose w lies on the axis of a. Then A = T and it follows from (|3.10|) that 
for all u £ V: 

T 3 (u) = T 2 ((u, v)w — (it, w)v) = T((u, v)v — (u, w)\v \ 2 w) 

= (it, v)t 2 w — (u, w)t 2 v — t 2 T(u), (3-17) 

where r = \v\. Thus (|3.15|l is satished with A = r 2 , so a is preharmonic with: 

C = 2F-t 2 , 

by (|3.16[> . Furthermore by (13.71) : 

AF = 2(t 2 - (n + l)F). 

Harmonicity equation (|I.9[) is therefore polynomial in the non-constant smooth 
function F (indeed, no non-trivial Killing field on a manifold of strictly negative 
Ricci curvature has constant length), and gathering coefficients yields: 

(n + 1 — p)q = 0, 
(n + 1 + (p - 2)r 2 )q = l-n, 
1 - n + 2qr 2 = 0. 

It follows (from the second or third equation) that q ^ 0, therefore (from the 
first) p = n + 1. Adding the second and third equations then yields 1 + r 2 = 0, 
contradicting the fact that v is spacelike. □ 

Proposition 3.5. Let a be an infinitesimal rotation of H n . Then a is preharmonic 
if and only if a is balanced. 

Proof. We may orthogonally decompose V = Lffi T w H n where L C V is the line 
through w. Suppose w lies on the axis of a. Then A — R, hence L C ker(A) and 
the restriction of A to T w H n is a skew-symmetric transformation of a Euclidean 
vector space. The argument of Proposition 13.21 then applies. □ 

We now show that there are no harmonic Killing fields of parabolic type. 

Proposition 3.6. Let a be a Killing field on H n . If a is harmonic then a is an 
infinitesimal rotation. 

Proof. Suppose a is harmonic. Since v £ im(A) and a is preharmonic it follows 
from (|3~15)) that A 2 (v) = \v. Now: 

RT(v) = R{t 2 w) = 0, 

and by (j3~10| and pTf) : 

TR(v) = TA(v) - T 2 {v) = (A(v),v)w - (A(v),w)v - t 2 v = 0. 

Therefore: 

A 2 (v) = R 2 (v) + T 2 (v) = R 2 {v) + t 2 v. 
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Since R may be regarded (by restriction) as a skew-symmetric cndomorphism of 
the Euclidean vector space T w H n , there is an orthogonal decomposition v = v + v 
where v E im(i?) and v E ker(_R). Then: 

A 2 (v) = R 2 (v) +t 2 v + T 2 ii. 

Therefore, since R 2 (v) E im(R), v lies in the A-eigenspace of A 2 if and only if: 

R 2 {v) = (\-t 2 )v and (A - t 2 )v = 0. (3.18) 

Case 1. Suppose -0 = 0. Then for all u E im(R): 

A(u) = R(u) + T(u) = R(u) + (u, v)w - (u, w)v = R(u). 

Since R restricts to an isomorphism of im(_R), this also shows that im(i?) C im(A). 
Proposition 13.41 ensures R ^ 0, so it follows from (|3.15l) that A is a non-zero eigen- 
value of R 2 . Then A < 0, since R is Euclidean skew-symmetric, hence v = by 
(|3.18|) . Therefore v — 0. Thus a is an infinitesimal rotation, with w on the axis. 

Case 2. Suppose i ^ 0. Since R 2 restricts to an isomorphism of im(i?) it follows 
from p,18[) that A — t 2 ^ 0, hence v = 0. So v E im(i?) is a non-null eigenvector 
of R 2 . Therefore v — V\e\ + • • • + V2 r &2r with respect to the orthonormal basis of 
(|3.8[) . We now introduce the index subset: 

/ = {< G {1, . . . ,r} : (v 2l -if + {v 2i ) 2 + 0}. 

It follows that all the u>i with i E I are equal; say a;, = p. Thus R 2 (v) = —p 2 v, and 
comparison with ([3.18J1 yields: 

\ 2 2 

A = r — p . 

Let J be the complementary index subset. Thus j G J if and only if i>2j-i = v 2 j = 0, 
and it follows that: 

T(e 2j -i) = = T(e 2j ). 

Hence: 

A(e 2 j-i) = R(e 2j -i) = Uje 2j , A(e 2j ) = R(e 2j ) = -u^j-i. 

So e 2 j-i 1 e 2 j E im(^4) and: 

A 2 (e 2j -i) = -uj 2 e 2 j-i, A 2 (e 2j ) = -u 2 e 2j . 

Therefore lj 2 = —A = p 2 — t 2 for all j E J, by (|3 . 1 5[) . In particular, p > r. Since 
R 2 has precisely two non-zero eigenvalues there is an orthogonal splitting: 

im(i?) = V © W, 

where V (resp. W) is the eigenspace for eigenvalue —p 2 (resp. t 2 — p 2 ); indeed, V 
(resp. W) is spanned by {e 2 i-±, e 2 i : i E 1} (resp. {e 2 j-i,e 2 j : j E J}), and v G V. 
Hi £ I then: 

T(e 2 i-i) = «2i-iw, 7 1 (e 2l ) = v 2i w, 

hence: 

A(e 2 i-i) = /fe 2l + u 2l -iw, A(e 2 ,) = -pe 2i _i + v 2l w. 

Therefore: 

A(v 2i e 2i -i - v 2l -ie 2i ) = p(v 2i -ie 2i -i + v 2l e 2l ), 
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so the vector <n = V2i-ie 2 i-i + vnen € V also belongs to im(A). Now: 

A 2 (a l ) = A(v 2l -i{pe 2l + v 2l -\w) + v 2 i{~pe 2l -i + v 2l w)) = ~p 2 ai + \ai\ 2 v. 

Since a, is an eigenvector of A 2 with eigenvalue A = t 2 — p 2 , we must have v = a,. 
We conclude that is 2-dimensional. In particular, this implies: 

ruj 2 - t 2 = p 2 + (r - l)(p 2 - r 2 ) - r 2 = r(p 2 - t 2 ) > 0. 

Therefore a is again an infinitesimal rotation, but with w no longer on the axis. □ 

The spherical and hyperbolic cases can henceforward be treated simultaneously. 
For fixed r with 2r < n + 1 we introduce the following twist equation: 

2cW + e(2nk - c)u 2 + 1 - n = 0, (3.19) 

where c = n + 1 — 2r and k = r — 1. If r = 1 then (|3.19|) has solutions only when 
e = — 1, with unique positive root ljo = 1. If r > 1 the discriminant A of (|3 . 1 9|) . 
viewed as quadratic in w 2 , is strictly positive, and the signs of the coefficients show 
that there is precisely one positive root uiq, which we refer to as the optimal twist. 
We claim that if e = 1 (resp. e = —1) then ujq < 1 (resp. ojq > 1). For, from the 
following upper bound: 

A = (2nk - cf + 8ck(n - 1) = (2nk + cf - 8ck < (2nk + c) 2 , 

we obtain for e = 1: 



UJ, 



, c - 2nk + 2nk + c 2c 1 . 

? < 171 = 731 = 777-^7- (3-20) 



4cfc 4cfc 2(r - 1 

So in the spherical case there is in fact the sharper estimate ujq < \j\[2. On the 
other hand, since n > c there is the following lower bound: 

A > (2ck - c) 2 + 8c 2 k = c 2 (2k + l) 2 , 

from which if e = — 1 we obtain: 



w a > = —^--^ = 1 + (3.21) 



2 ^ 2nfc-c + c(2fc+ 1) 

~\~ck 2c ' 2c 

It was noted (Corollary 13. 3[) that on odd-dimensional spheres the harmonic 
Killing fields of maximal rank are constant multiples of Hopf fields, which are 
precisely the spherical Killing fields of constant length. By contrast, there are no 
non-trivial hyperbolic Killing fields of constant length. The following result deals 
with all other cases; ie. Killing fields of non-constant length. Recall the generalised 
Hopf field S r defined in fl3U). 

Theorem 3.7. Let a be a Killing field of non- constant length on the non-flat space 
form M . Then a is harmonic if and only if a is congruent to cJoS r for some 
< r < (n+ l)/2, with r > 1 if M — S n . Furthermore a is metrically unique, with 
metric parameters p = n + 1 and: 



n)/2, ifr = \ 

if r > 1 



o 2 . . (t) 
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Proof. If a is harmonic then it follows from Propositions 13.11 13. 2\ 13.51 and 13.61 that 
in the spherical case a is balanced, and in the hyperbolic case a is a balanced 
infinitesimal rotation. Then equation (|3.15l) holds with A = —to 2 , where oj is the 
twist. By (|3.16[) the spinnaker of a is: 

C = uj 2 - 2eF. 

Furthermore by Q3.7P : 

AF = 2{e{n + l)F -ruj 2 ). 

Harmonicity equation (|1.9p is therefore polynomial in the non-constant function F, 
and the coefficients yield: 

(n + l-p)q = 0, (3.22) 
(e(n + 1) + - 2r)uj 2 )q = e(l - n), (3.23) 
2((r - l)p + rq)u 2 =e(n-l). (3.24) 

It follows from (f3T23|) that q ^ 0, whereupon p = n + 1 by (|3^2| . Adding (|3T23|) 
and (|3.24j) then yields: 

2(r-l)w 2 + g(cj 2 +e) = 0. (3.25) 

If r = 1 then s = -1 and lu = 1, and it follows from (13T2"4]) that q = (1 - n)/2. If 
r > 1 then (13.25)) rearranges to ([[]) for g, substitution of which into (I3.23[) yields 
the twist equation p. 191) for uj. □ 

It follows from ([]]) that all harmonic Killing fields a of non-constant length have 
metric parameter q < 0. (Recall that in the hyperbolic case lvq > 1 when r > 1.) 
In the spherical case, plugging the upper bound (|3.20|) for ljq into ([[]) yields the 
following lower bound for q: 

2k 2k 2r-2 . , 

9 = "ITIM? > _ TT2fc = -27^1' (3 - 26) 

Thus — 1 < g < 0, in contrast to harmonic conformal gradient fields, where q ^ — 1. 
Since w 2 < 1 it follows that go; 2 > —1, which since uiq = HcHcc shows that a 
is g-Riemannian (|1.5|) . In fact, when (|3.26p is combined with p.20p the following 
sharper estimate is obtained: 

2 1 1 1 

quj n > > > — , since r ^ 2. 

y 2r- 2 2r- 1 3' 

So cr(M) sits comfortably inside the g-Riemannian ball bundle. On the other hand, 
bearing in mind that n > 3 if a has non-constant length: 

n n + 1 

It therefore follows from Proposition 11.11 that w 2 > 1/n. We conclude that: 
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These bounds are tightest when n is even and a has maximal rotational rank 
r = n/2. For example, applying (|3.26j) and (|3.27|) to the unique (up to congruence) 
non-trivial harmonic Killing field on S yields: 

-1/3 <?<(), 1/4<w 2 <1/2, 

in comparison to the precise values: 

q = (V73 - 13)/8, w 2 = (V73 - 7)/4. 

In the hyperbolic case, if n ^ 4 and a has rotational rank r > 1 then since ujq > 1 
it follows from ((fj that q < 2(1 - r). In particular, g < -2. From (|3.21[) : 

w 2 1 n + c 

ujq — 1 1-1 /u; 2 rz — c ' 

which by ([f]) yields a lower bound on q: 

2(5(1 - r) <g< 2(1- r), where (5 = — — > 1. (3.28) 

n — c 

These bounds are again tightest when n is even and r = n/2. For example, applying 
(13.281) and (13.21)) to the unique (up to congruence) harmonic Killing field on H 4 
with r = 2 yields: 

-10/3 < q < -2, cj 2 > 5/2, 
in comparison to the precise values: 

g=-(V73 + 13)/8, w 2 = (V73 + 7)/4. 

4. Loxodromic Fields 

Suppose a, € V and let a, /? be the restrictions to M of the corresponding 
metrically dual linear functionals. First define for all x G M; 

K(x) =a(x)b- /3(x)a. (4.1) 

Then K is the restriction to M of a skew-symmetric transformation of V, hence a 
Killing field. If A, _B are the conformal gradient fields with poles a, b then by (|2 . 3|) : 

K{x)=aB - PA. (4.2) 

We refer to if as the eiementary Killing Held determined by the pair (a, b). From 
(|2.6[) . the covariant derivative of K is particularly simple: 

V X K={A,X)B-{B,X)A. (4.3) 

It is also useful to note from (I2.3[) the following generalisation of (|2.4I) : 

(A,B) = {a,b)-eap. (4.4) 

Now let ^ = {ai, &i, . . . , a r , b r } be a spacelike orthonormal subset of V. Let Ki 
be the elementary Killing field determined by (a^, and define: 
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From Section [3l every Killing field on S n , or infinitesimal rotation on H n , of ro- 
tational rank r has such an expansion. A vector field a on M is then said to be 
loxodromic if: 

a = R + C, 

where C is a non-trivial conformal gradient field whose pole c is orthogonal to 8£. 
We say that a is properly loxodromic if R is balanced and n — 2r; in the hyperbolic 
case this implies that c is timelike. Then R and C have precisely the same zeros, 
and the angle between u{x) and C{x) is constant. If n = 2 then every loxodromic 
field is properly loxodromic. 

Let 7 be the restriction to M of the linear functional metrically dual to c. It 
follows from (|4~2]l and (|4~4]) that: 

(Ki, C) = ai {Bi, C) - /3i{Ai,C) = -ea^j + e^a. a = 0. 

Furthermore since is spacelike: 

(K i ,K j ) = 8 tj ia? + tf), (4.5) 

where 5ij is the Kronecker symbol. So if a is loxodromic then (summing over i): 

M 2 = Lo?\Ktf + |Cf = tof{af + tf) + M - e 7 2 , (4.6) 

where (i = (c, c). Note that a vanishes precisely on the axis of R; ie. the totally 
geodesic submanifold M n 3% C M of codimension 2r. 

Lemma 4.1. Let a be a loxodromic vector field on a non-flat space form. If a is 
p-preharmonic for some p then a is properly loxodromic. 

Proof. From and (|2~%|) : 

V*Va = e{n-l)R + eC. (4.7) 

Therefore <r is 0-preharmonic if and only if n = 2, in which case a is properly 
loxodromic. Now suppose p ^ 0. Introducing the gradient vector fields: 

d = diAi + PiBi, 

it follows from (|4.6p that (summing over i): 

VF = u?Gi - e7C. (4.8) 

From 

(G i ,K j )=8 ij (a j 0i-ai0 j ) = O. 
Therefore VF is orthogonal to i?. From (|4.3[) and (|2.6[) (summing over i): 

V vf ct = uj i {A i ,VF)B i - Ui {Bi, VF)A l - ejWF. 
Then by (I4.8[) and (|4.4[) (summing over all repeated indices): 

V VF cr -I- sjVF = Uibjf {otj{5ij - Eauatj) - £aiPf)Bi + us l a i ^ 2 B l 

- uiiU^PjiSij - efiipj) - ef5ia])Ai - utfi-fAi 
= e(/i + ew? - M^Wj-K,-, by gU). (4.9) 
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Since VF is orthogonal to R, it follows from ()4.7j) . (|4.9|) and (|1.4|) that if a is p- 
preharmonic then R is balanced and V-F is pointwise collinear with C . The latter 
is equivalent to the pointwise collinearity of G = Gi and C. Using (14.4)) we first 
note that for all i, j: 

(Gi,Gj) = (a 2 + - s(a 2 + /?/)), 

and then calculate (summing over repeated indices): 
|CAG| 2 = \C\ 2 \G\ 2 - (C,G) 2 

= (a 2 + tf) {n - e 7 2 - BAt (a? + /?/)) 

Since -£Q vanishes only on a set of measure zero, if e = — 1 then C AG vanishes only 
if fi < 0; ie. c is timelike. So there exist vectors d±, . . . , d n ~2r € V, which in the 
hyperbolic case are spacelike, such that {a\, b\, . . . , a r , b r , d\, . . . , d n ^2r, c/ ^/sjl} is 
an orthonormal basis of V. Hence, if Si : M R is the restriction of the covector 
metrically dual to di then there is the quadratic relation: 

E ■=!(«/ + P?) + S 2 + --- + S 2 _ 2r + 7 2 ^ = e, (4.10) 

which rearranges to: 

\C\ 2 -e^J2 3 \K,\ 2 =e^S 2 + --- + 5 2 _ 2r ). 

Therefore if C A G = then n = 2r. □ 

Proposition 4.2. Let a be a loxodromic vector field on the non-flat space form 
M n . Then a is p-preharmonic for some p if and only if n = 2. In this case a is an 
eigenfunction of the rough Laplacian, and a is preharmonic with spinnaker: 

Proof. Suppose that a is properly loxodromic. Since n — 2r, relation ()4.10|) reduces 
to (summing over i): 

M (a 2 + /? 2 )+ 7 2 = £ M , (4.11) 
which when differentiated yields: 

MG + 7 C* = 0. (4.12) 
Furthermore since R is balanced (|4.9p simplifies: 

V VF a = e{fi + eu 2 - \a\ 2 )R - e 7 VF. 
Application of (|4TTT|) to jUB]) yields: 

/i|a| 2 = ( M + eu 2 )(»-s 1 2 ). (4.13) 
Also application of (|4.12|) to f|4T8]) yields: 

fiVF = -e(fi + ew 2 ) 7 C, (4.14) 

and then by ([4TT3| : 

M7 VF = -e(/x + eu; 2 ) 7 2 C* = + ew 2 - |er| 2 )C. 
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This yields the further simplification: 

V VF (T = e(fi + ecu 2 - |cr| 2 )cr. 

It then follows from (|4.7j) that a is p-preharmonic if and only if n 
is now a consequence of Lemma 14.11 

Define vector fields cto and <j\ on H 2 as follows: 

a (xi,x 2 ,x 3 ) = (-x 2 ,xi,0), 
<ji(x 1 ,x 2 ,x 3 ) = (-xix 3 ,-x 2 x 3l l -x 2 ). 

Then cto (resp. o~\) is an infinitesimal rotation (resp. conformal gradient) with axis 
(resp. pole) (0, 0, 1), and any loxodromic field on H 2 is congruent to ua a + ^—[ia\. 
By Theorems 12.21 and 13.71 up to congruence uo (resp. <j\) is the unique harmonic 
Killing (resp. conformal gradient) field on H 2 , and both fields are metrically unique, 
with metric parameters (p,q) = (3,-1/2). Furthermore, when H 2 is regarded as a 
Riemann surface Co and <j\ are conjugate. The associate family: 

Jz? = {sin t do + cos t cti :tel} 

is therefore also harmonic, and loxodromic (with the exception of do and a%, which 
our definition excludes). 

Theorem 4.3. Let a be a loxodromic vector field on the non-flat space form M . 
Then a is harmonic if and only if M = H 2 and a is congruent to an element of 
Jzf '. Furthermore a is metrically unique, with metric parameters (p,q) — (3,-1/2). 

Proof. In Proposition 16. II we will show that when n = 2 the harmonicity equation 
(|1.3j) for a conformal field simplifies: 

e(l + 2F)(1 + 2qF) + 2q((p ~ 2)F - 1)C = 0, 

where C is the spinnaker. It follows from Proposition 14.21 that this equation is 
polynomial in F, with leading coefficient (3 — p)q. Therefore p = 3, since there is 
clearly no solution if q = 0. The terms of lower degree then yield: 

2q{n + elj 2 ) = 1 = -({! + ew 2 ), 

from which q = —1/2 and e = — 1. Then \i < and ui 2 — \i = 1. □ 

5. Conformal Dipole Deformations 

First let M = S n . Pick any w G S n and a G T w S n . The Pontryagin vector 
Geld [TT] determined by this point and tangent vector is constructed by parallel 
translating a along the geodesies from w, to obtain a vector field of constant length 
with a singularity at —w. The singularity may be smoothed off to a zero by a 
natural scaling factor — the square of the cosine of half the geodesic distance from 
a — resulting in a smooth vector field 5 with a point dipole at —w. Analysis of this 
procedure, which we omit for brevity, produces the following expression: 

8(x) = i((l + i/)(x))a - a(x)(x + w)), for all x e S n , 



= 2. The result 
□ 
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where ip, a are the metric duals of w, a. It follows from (|4.2j) and ()2.3|) that: 

25 = A — T, 

where A is the conformal gradient with pole a, and T is the elementary Killing field 
determined by the pair (a,w). If W is the conformal gradient with pole w then: 

T = aW-il)A. (5.1) 

The dipole may be repositioned from — w to w by reversing the sign of a, and any 
constant rescaling of the field may also be incorporated into a. We therefore define, 
for M = S n or M = H n , the dipole field determined by w G M and a G T W M to 
be: 

<j = T + eA. 

The geodesic through w tangent to a is the dipole axis. In the hyperbolic case, A is 
a conformal gradient field without zeros, and comparison of (I5.1[) and (|3.10[) shows 
that T is an infinitesimal translation whose axis coincides with that of the dipole 
and meets the equator of A orthogonally at w. 

More generally, suppose a € T W M is a unit vector and define for r, r e R: 

<T = rT + rA. (5.2) 

We refer to a as a dipole deformation Geld. Bearing in mind that (w, w) = e and 
(a, to) = it follows from (|4.4[) : 

H 2 =r 2 |T| 2 + 2rr(T,A) + r 2 |A| 2 

= (rV>-r) 2 +£(r 2 -r 2 )a 2 . (5.3) 

If r = r and M = S n , or r = — r and M = H n , then er is a dipole field and it 
follows from (|5.3p that cr has a single zero, at u>, as expected. Otherwise, depending 
on the relative magnitudes of r and r, the point dipole at w splits into a pair of 
zeros on the dipole axis, with midpoint w; or the zero set of a is a codimension- 
1 submanifold of the equator of A (ie. codimension-2 in M) with two connected 
components, placed symmetrically with respect to w. 

Lemma 5.1. Let a be a dipole deformation field on the non-flat space form M n , 
that is neither a conformal gradient nor a Killing field. Then a is p-preharmonic 
for some p if and only if n = 2. In this case a is an eigenfunction of the rough 
Laplacian, and a is preharmonic with spinnaker: 

( = e(r 2 +t 2 - 2rrV> - M 2 )- 

Proof. By assumption t^O and r ^ 0. From (|5.3p : 

VF = t{ttP -r)W + e(r 2 - r 2 )aA. (5.4) 

Therefore by (|^6]l . ([44]) and fBTT]) : 

V VF fi + eraVF = r(A, VF)VK - t(W, VF)A 

= sr 2 (r - Ttp)(A + tjiT) + t(t 2 - r 2 )a(£l^ - aT), 
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which after further extensive calculation yields: 

V VF cr = (r 2 - r 2 )a 2 a + er 2 (l - ^ 2 )a. (5.5) 
Now from (H) and (j3~T]) : 

V*V<r = e(n - 1)tT + rA. (5.6) 

Hence a is p-preharmonic precisely when n = 2. The expression for £ follows by 
comparison of (|5.5[) with (|5.3I) . □ 

Theorem 5.2. Lei a be a non-trivial dipole deformation field on a non-flat space 
form. Then a is harmonic if and only if a is a harmonic Killing field or a harmonic 
conformal gradient field. 

Proof. By Lemma I5TT1 it suffices to consider the case n — 2. Anticipating Proposi- 
tion [6T] the harmonicity equation then simplifies: 

e(l + 2F)(1 + 2qF) + 2q((p - 2)F - 1)C = 0, 

where ( is the spinnaker. It follows from (|5.3p and Lemma 15.11 that this is a 
polynomial P(a, ip) — 0. The highest order (quartic) terms are: 

P 4 (a, if,) = (3 - p)q{e{T 2 - r 2 ) 2 a 4 + 2r 2 (r 2 - r 2 )a 2 ^ 2 + er 4 ^ 4 )- 

Since q ^ and a is non-trivial, p — 3. The cubic terms then reduce to: 

P 3 {a,ip) = 6qrr((r 2 - r 2 )a 2 ip - eT 2 iP 3 ) , 

the vanishing of which forces r = 0orr = 0. □ 

6. Conformal Vector Fields in Dimension Two 

Suppose that a is a conformal vector field on M = S 2 or M — H 2 . Let J 
be one of the two Kahler structures on (M,g) (the choice is immaterial). The 
conformal transformations of M are then precisely its biholomorphic mappings, so 
a is a holomorphic vector field. Furthermore, a may be decomposed as the sum of 
a Killing field and a conformal gradient field, so it follows from (|2.8p and (|3.ip that 
a is an eigenfunction of the rough Laplacian, with eigenvalue e. Our first result 
uses these observations to recast the harmonicity equation (jl.3|) for a. 

Proposition 6.1. Let a be a conformal vector field on M . Then a is preharmonic; 
and a is (p, q) -harmonic if and only if: 

e(l + 2F)(1 + 2qF) + 2q((p - 2)F - 1)C = 0, 
where £ is the spinnaker of a. In particular, if a is (p, q) -harmonic then q =/= 0. 
Proof. Since a is a holomorphic vector field we have: 

V .( J<J ) = + K Jo~] - V Jct (t + J[a, a] = V Ja a, 

because L a J = 0. The components of Vc may therefore be obtained: 
(V a a,a) = (VF,a) - (V a (Ja),Ja) = (V Ja a,Ja), 
(VF,Ja) = {o-,V Ja a) = (a,\(Ja)) = -{\a, J a). 
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Now: 

2FVF = (VF,a)a + {VF,Ja)Ja. 

Hence: 

2F(V vf( t, Ja) = (VF, a)(V a a, Ja) + (VF, Ja)(V Ja a, Ja) = 0. 

Because the zeros of Ja are discrete, this implies that V VF <7 is pointwise collinear 
with a. Therefore, since also a is an eigenfunction of the rough Laplacian, a is 
preharmonic. We now observe: 

2F|VF| 2 = (VF,cj) 2 + (VF, Jcr) 2 

= °? + ^ J °f + & Ja a, a) 2 + (V Ja a, Ja) 2 ) 

= 2F 2 \Va\ 2 . 

By comparison with (|1.7j) the spinnaker of a must therefore be: 

C=^|Va| 2 . 

It then follows from the Weitzenbock identity f|l .8[) that: 

AF = £ | ( r| 2 -|Va| 2 -2(eF-C), (6.1) 

which yields the stated simplification of harmonicity equation (|1.9j) . The equation 
clearly cannot be satisfied if q = 0. □ 

We will now represent: 

a = K + C, 

where K is a Killing field and C a conformal gradient. To do so, first choose w £ M 
and an orthonormal basis (a, b) of T W M. Let a,f3,ip denote the restrictions to M 
of the linear forms on V metrically dual to a, 6, w respectively, and define Killing 
vector fields on M: 

R = aB — j3A, T = aW — ipA, 

where A,B,W are the conformal gradients with poles a,b,w respectively. In the 
hyperbolic case R is an infinitesimal rotation about w, and T is an infinitesimal 
translation through w, therefore: 

K = ujR + tT, 

for some uj, t ^ 0. In the spherical case the situation is simpler, for by choosing w 
to lie on the axis of K it may be assumed that t = 0. Now let c be the pole of C, 
and let 7 be its metrically dual covector. We locate c cylindrically: 

c = rsa + rtb + hw, 

where h £ R, r ^ and s 2 + 1 2 = 1. Finally we note the quadratic relation: 

Q(a,/3,V) = a 2 +(3 2 + evjj 2 -e = 0, (6.2) 

which expresses the equation of M with respect to the "coordinates" a,/3,ip. 

Lemma 6.2. If a is a conformal vector field on M 2 then: 

|<r| 2 = t 2 +r 2 + e(uj 2 + h 2 ) + 2(urt + erh)a - 2rs(w(3 + rip) - e(ujip - er/3) 2 ~ ej 2 . 
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Proof. We have: 

H 2 = oj 2 \R\ 2 + t 2 \T\ 2 + |Cf + 2ut(R, T) + 2u(R, C) + 2r(T, C). 
Each inner product of conformal gradients may be computed from (|4.4|l : 
\R\ 2 = a 2 +/3 2 , \T\ 2 = ea 2 + V> 2 , \C\ 2 = r 2 + eh 2 , 
(R,T)=/3ip, (R,C) =r(ta- s/3), (T, C) = eha — rsip, 
and the result follows after applying the relation (|6.2p . □ 
Lemma 6.3. If a is a conformal field on M 2 then the spinnaker of a is: 

C = {Lu^-erp) 2 + 7 2 . 

Proof. From Lemma 16.21 

VF = {tort + erh)A - 2rs{uoB + tW) - e{u)ip - st/3)(ujW - erB) - e^C. 
Therefore by flU) and (f2~T0j) : 

AF = 2e(wri + erh)a - 2ers(ujf3 + Tip) - 2{ujip - erfi) 2 - 2j 2 
+ e\u:W -erB\ 2 +s\C\ 2 
= e\a\ 2 - 2{u)ip - erf3) 2 - 2 7 2 , 
by (|4.4[) and comparison with Lemma 16.21 The result follows from (|6.1[) . □ 

It follows from Proposition 16.11 and Lemmas 16.21 and 16.31 that a is harmonic if 
and only if: 

P(a,/3,TP) = 0, 

where P is a quartic polynomial. This is equivalent to the algebraic problem of P 
vanishing modulo the polynomial Q of (16. 2|) : thus if Pk denotes the terms of P that 
are homogeneous of degree k we require Pk = {QS)k for all ^ k ^ 4, for some 
(necessarily quadratic) polynomial S. Now: 

P 4 (a,/3,V) =eq(3-p)( 2 1 P 3 (a,/3,V>) = ff(p- 4)r?C, 
where ?y denotes the linear part of |cr| 2 , and g ^ 0. On the other hand: 

(Q5) 4 (o, /8, 1>) = (a 2 +P 2 + £^ 2 )S 2 (a, (3, i/>), (6.3) 
and because Q\ = 0: 

(QS) 3 (a, [3, V) - (a 2 + (3 2 + £ip 2 )S x {a, f3, ip). (6.4) 

We consider the spherical and hyperbolic cases separately. 

Theorem 6.4. No non-trivial conformal field on the 2-sphere is harmonic. 

Proof. In addition to the simplification r = 0, by rotating the basis (a, b) in T W S 2 
if necessary it may be assumed that s = and t — 1. Then: 

C = u 2 ip 2 + 7 2 = (oj 2 + h 2 )^ 2 + 2rhf3ip + r 2 f3 2 . 

In view of Theorem 12.21 we assume w^O. Furthermore, a is loxodromic (resp. an 
element of a dipole deformation) when r = (resp. h = 0), so by Theorem 14.31 
(resp. Theorem 15. 2p we assume r^O and h ^ 0. Since £ is independent of a, no 
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match of P4 with (|6.3[) is possible unless P4 vanishes identically; hence p — 3. Now 
rj = 2ura, therefore P3 has a non-vanishing afiip term, which means that no match 
with (|6.4p is possible unless P3 vanishes identically, which is also impossible. □ 

Theorem 6.5. If a is a non-trivial harmonic conformal field on the hyperbolic 
plane then a is congruent to an element of the associate family Jzf . 

Proof. If lu = h = then a belongs to a dipole deformation; so by Theorem l5.2l we 
assume w 7^ or h ^ 0. The coefficient of ip 2 in £ is lS 2 + h 2 , so S2 contains the 
term — (a; 2 + h 2 ) 2 ip 2 . The coefficient of a 2 in £ is r 2 t 2 , so S2 also contains the term 
r A t A a 2 . Therefore (QS)a contains the term: 



whose coefficient is strictly negative. However the coefficient of a 2 ip 2 in P4 is 
4r 2 i 2 ft 2 . It follows that P 4 vanishes identically, hence p — 3. Now, inspecting the 
coefficients of a 3 ,/3 3 and tp 3 in P 3 yields: 

Si = 2q((rh - ujrt)r 2 s 2 a + wrs(r 2 + r 2 i 2 )/3 - rrs(w 2 + h 2 )^) . 

However the coefficients of aip 2 and (3ip 2 in P3 are, respectively: 

2q(rh-ujrt)(uj 2 + h 2 ), 2qujrs(uj 2 + h 2 ), 

and comparison with the corresponding coefficients in (QS)3 yields equations: 

(rh - Lort)(u? + h 2 + r 2 s 2 ) = = Lors{uj 2 + r 2 + r 2 t 2 + ft 2 ), 

which reduce to: 

t/i — Lurt = = urs. 

If u> = then ft ^ 0, so t = 0, and c is therefore a conformal gradient. If r = 
then either ft = in which case a is a Killing field, or r = in which case a is 
loxodromic. Finally, if s = then: 



2F = t 2 - uj 2 + r 2 - h 2 + C = + C, say. 
After some gathering of terms: 

P 2 (a,/3,VO = -(l + (0 + %)C 

Since £ is now independent of a, P2 is not a multiple of Q2 and therefore vanishes 
identically: 



(r 4 < 4 + ( W 2 + ft 2 ) 2 )aV, 




(6.5) 



with: 



(0 + 3)q = -l. 

Now P has no linear terms, so it remains to collate the constants: 

P {a,f3^) = -{l + 6){l+qd). 



(6.6) 
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The only solution of Po = consistent with (|6.6|) is: 

e = -i. (6.7) 

If a is an infinitesimal translation or of parabolic type (ie. r 2 ^ uj 2 ) then it follows 
from (|6.5|) that 9^0, contradicting (|6.7[) . If cr is an infinitesimal rotation then 
choosing for convenience w to be the axis, so that r = 0, it follows from (I6.5[) that 
r = 0; thus cr is again loxodromic. The result now follows from Theorem 14.31 □ 

7. Quadratic Gradient Fields on Spheres 

Let M — S n with n > 1. Let Q: V — ► V be a symmetric linear transformation, 
and let £ : M — »• R be the restriction of the associated quadratic form: 

£(a;) = (Q(x),x), for all £ € M. 

Now define: 

<7=±V£. 

If (oi, . . • , a n +i) is an orthonormal Q-eigenbasis of V: 

Q{ai) = Xidi, Xi € R, 
and cti : 5™ — >• R are the restrictions of the metrically dual covectors (|2.1[) , then: 

e = E l A I a 2 . 

Therefore: 

<T = E i Aia i i4 j) (7.1) 

where Ai is the conformal gradient field with pole <Zj . Applying (|2.3|) to (IT. ip yields 
the following coordinate-free representation: 

£7 ( a; ) = Kai(x)(ai ~ a>i(x)x) = Q(x) ~ £(x)x. 

In particular, this shows that the zeros of a are precisely the (unit) eigenvectors 
of Q. Thus, a has at least 2n + 2 zeros, and these are isolated if and only if 
the eigenvalues of Q are simple. It should also be noted that Q is not uniquely 
associated to a, because of the relation: 

a i H 1" = !• 

In particular, the eigenvalues Ai do not determine a up to congruence. The confor- 
mal representation (|7.1[) is also non-unique, because of the differentiated relation: 

ai Ax + ■ ■ ■ + On+xAn+x = 0. (7.2) 

However, (|7.1I) does provide computational assistance, using results of Section [5] 
We first compute the divergence of a (summing over i): 

div(cr) = X l (a l divA l + \A t \ 2 ), by fl^TO) 

= A i (l-(n + l)a i 2 ), bv (I2T71) and (l2~l) 

= tr(Q)-(n + l)£. (7.3) 

Also, using the general Riemannian identity: 

V*V(/A) - fV*VX - 2V V ,A + (Af)X, 
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we compute using (|2.8|) . (|2.6[) and (|2.7p (summing over i): 

V*Vcr = \ i ot i A l + 2\. l a l A i + nX^A., = (n + 3)a. (7.4) 



Thus a is an eigenfunction of the rough Laplacian. It follows from (|7.4[) and (jl.4|) 
that when p = the harmonicity equation (|1.3j) reads: 



gAF + n + 3 = 0, 

which contradicts the Divergence Theorem. Therefore no non-trivial harmonic a 
has metric parameter p = 0. Hence if er is harmonic then a is preharmonic. 

Now for any positive integer to let £ m : M — >• M be the quadratic form associated 
to Q m = Q o ■ ■ ■ o Q (m iterations) : 

U = E,(^) m a, 2 - 
and let a m be the corresponding gradient field: 

a m = iV£ m = ^(A*)" 1 ^, (7.5) 
with coordinate-free expression: 

o m {x) =Q m {x)-i m {x)x. (7.6) 
It follows from (|4.4[) that for all fc,TO € N (summing over i,j): 
(<Jk,v m ) = (X l ) k (X j y n a l a J (A l ,A J ) 

= (\i) k+m af - (X^a^X^af = £ k+m - (7.7) 

In particular: 

H 2 = 6-C 2 - (7.8) 

Then: 

VF = cr 2 - 2£er. (7.9) 

We note that in addition to vanishing at the zeros of a, VF also vanishes at the 
midpoints ay = (a; + dj )/v / 2, and by (EH]): 

kKOI 2 = |(A 2 + A 2 ) - i(Ai + A,) 2 = i(A, - A,) 2 . 

Therefore: 

||cr||oo = maxi|Ai - Xj\. (7.10) 
Of course, if all eigenvalues of Q are equal then a is trivial. 

Proposition 7.1. Let a be a non-trivial quadratic gradient field on S n . Then a is 
preharmonic if and only if Q has precisely two distinct eigenvalues. 



Proof. It follows from (|7.9p . (|7.ip and (12 .6[) that (summing over i,j): 

7 VF l 



V VF a = X i (VF,A i )A i - X t af VF 



= Xi XjCij{Ai,Aj)Ai - 2;\\,A,<», .1, £VF 

which by P~4"|) : 

= ^3 - 3^ 2 + (4£ 2 - (7.11) 
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It therefore follows from (|7.11|) that a is preharmonic if and only if: 

<x 3 - 3^ 2 = fa, (7.12) 

for some smooth function / : M — > R. Since a is non-trivial there exist j, k such 
that Xj 7^ Afe. Set x = dj 6 M and X = a k € T X M . Covariant differentiation of 
equation (|7.12j) along X yields: 

V x cr 3 - 3Aj V x (t 2 = f(aj)W x a, 

bearing in mind that criflj) = 0. From the conformal representation (|7.5[) , using 
(f2~6]) and (|2~3]l we obtain for all m E N: 

V x <r m = ^(^((X^iKOMiK-) - = ((X k ) m - (X,) m )a k . 

Therefore if a is preharmonic: 

A, 3 - X? - 3A,(A 2 - X?) = (X k - Xj)f( aj ), 

and cancelling the non-zero common factor yields: 

f{a ] ) = Xl-2X J X k -2X]. 

Now if Xi is a third distinct eigenvalue then: 

f( aj ) = X} -2A,A,-2A 2 , 

and eliminating fiflj) yields: 

X k + Xe = 2Xj. 

By symmetry (ie. taking x — a k and X = dj): 

Xj + Xg = 2Afe, 

which forces the contradiction Xj = X k . 

Conversely, if Q has precisely two distinct eigenvalues then it follows from (17.1[) 
and (|7.2[) that a may be represented: 

a = X{ot\A\ H hMr), (7.13) 

for some r < n + 1, where A is the difference of the eigenvalues. The corresponding 
potential is: 

e = AK 2 + --- + ar 2 ). (7.14) 

Then for all m G N: 

Cm = A m -^, cr m = A m_1 er, (7.15) 
and (|7.11[) therefore reduces to: 

V VFt r = (A — 2£) 2 <7. (7.16) 
Thus ex is indeed preharmonic. □ 
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The quadratic gradient field a on S n represented by (17.13|) is congruent to AS,- 
where: 

S r (x) ((1 s r )xi, . . . , (1 s r )x r , S r X r ^i , . . . , S r X n +\) , 
for all a; = (x\, . . . , x n +i) G S' n , using the abbreviation: 

S r = X^ ~t~ ' ' ' ~\~ x r . 

The integral curves of a are great circle arcs from a great (n — r)-sphere to the 
orthogonal great (r — l)-sphere; and by (|7. 10[) the maximum speed is |A|/2, which 
is attained midway along each arc. For r 3 let Ao be the unique positive root of: 

(r - 2)A 4 + 2(r 2 - 5)A 2 - 8(r + 1) = 0. (7.17) 

If A is the discriminant of (|7.17l) then one obtains by elementary means: 

4(r 2 - 2) 2 < A < 4(r 2 - l) 2 , 

provided r ^ 4, with the upper bound in fact valid if r ^ 3. These yield bounds: 

7^2 < A " < t4t <" 8 ' 

Theorem 7.2. j4 non-trivial quadratic gradient field a on S" is harmonic if and 
only if n ^ 5 is odd and a is congruent to AoS r with r = (n + l)/2. Furthermore 
a is metrically unique, with p = r + 1 and: 

9 l + r + (A /2)2- W 

Proof. We noted earlier that if cr is harmonic then ct is preharmonic. Therefore by 
Proposition 17.11 it may be assumed that £ and a are of the form (|7.14[) and (|7.13[) 
respectively. It then follows from (|T. 15[) that (|7.8[) and (|7.9p simplify: 

|a| 2 = A£-£ 2 , VF = (A - 2^)tr. 

Therefore by ([2~T0)) and ([73]) : 

AF = - div VF = (2$ - A) div a + 4|cr| 2 

= -rA 2 + (n + 5 + 2r)A£ - 2(n + 3)£ 2 . 

Furthermore, from (|T. 16[) the spinnaker of a is: 

C = (A - 2£) 2 . 

Harmonicity equation (|1.9I) is therefore polynomial in £. From the quartic terms: 

(n + 3- 2p)q = 0. 

If q = then (jl.9|) reduces to: 

n + 3 + (2 - r)pA 2 + (n + 3 + 2(r - 3)p)A£ - (n + 3 + (n - 5)p)£ 2 = 0, 
from the constant and linear terms of which: 

(r - 2)pA 2 = n + 3 = 2(3 - r)p, 
forcing a disparity of sign. Therefore q ^ and 2p — n + 3. From the cubic terms: 

n + 1 - 2r = 0. 
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Finally, the constant and linear/quadratic terms, respectively, yield: 

({r + l){r-2)+rq)\ 2 = 2(r + l), (7.19) 
(4(r + l) + A 2 ) g = 2(r + l)(2-r). (7.20) 

Eliminating q yields equation (|7.17p for Ao, and (|7.20j) rearranges to Q. □ 

The metric parameter q of Theorem 17.21 is manifestly negative, with: 

2q>2-r. (7.21) 

The metric parameters (p, q) for a therefore satisfy: 

1 — r p 

Proposition 11.11 then implies that: 

k^ = \w\\l> 1 - 1 



p — 1 r 

This allows the following modification of (|7.18[) : 

which is in fact a more accurate estimate when r > 8. Combining (|7.2ip with the 
upper bound for Ao yields: 

2-r 1 



2(r-2) 2 

Therefore in all cases a is comfortably q-Riemannian. Furthermore, from (|7.19[) 
and HOP]) : 

4rq= (r-2)A 2 -2(r-l) 2 , 
which when combined with (|7.18[) yields tighter bounds for q: 

3-2(r-l) 2 4-2(r-l) 2 ,„ . 

— *r^ <q< — at^' (7 - 23) 

provided r ^ 4, the upper bound being valid if r ^ 3. (Use of (17.221) would 
marginally improve the lower bound when r > 8.) For comparison, the precise 
values of Aq and q for low dimensions may be tabulated: 



n 


r 


P 


Q 


A 2 /4 


5 


3 


4 




a/3-1 


7 


4 


5 


x/201 - 29 


V201- 11 


16 


8 


9 


5 


6 


VM- 13 


n/34- 5 


5 


3 
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